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ABSTRACT 

The equations governing general relativistic, spherically symmetric, hydrodynamic 
accretion of polytropic fluid onto black holes are solved in Schwarzschild metric to 
investigate some of the transonic properties of the flow. Only stationary solutions are 
discussed. For such accretion, it has been shown that real physical sonic points may 
form even for flow with 7 < | or 7 > |. Behaviour of some flow variables in 

the close vicinity of the event horizon are studied as a function of specific energy and 
polytropic index of the flow. 
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1 INTRODUCTION 

Investigation of accretion processes onto celestial objects 
was initiated by Hoyle & Littleton (1939) by computing the 
rate at which pressure-less matter would be captured by a 
moving star. Subsequently, theory of stationary, spherically 
symmetric and transonic hydrodynamic accretion of adia¬ 
batic fluid onto a gravitating astrophysical object at rest 
was formulated in a seminal paper by Bondi (1952) using 
purely Newtonian potential and by including the pressure 
effect of the accreting material. Later on, Michel (1972) 
discussed fully general relativistic polytropic accretion on 
to a Schwarzschild black hole by formulating the govern¬ 
ing equations for steady spherical flow of perfect fluid in 
Schwarzschild metric. Following Michel’s relativistic gener¬ 
alization of Bondi’s treatment, Begelman (1978) discussed 
some aspects of the critical points of the flow for such an 
accretion. Spherical accretion and wind in general relativ¬ 
ity have also been considered using equations of state other 
than the polytropic one and by incorporating various ra¬ 
diative processes (Shapiro, 1973a,b, Blumenthal & Mathews 


1976, Brinkmann 1980). Recently Malec (1999) provided the 
solution for general relativistic spherical accretion with and 
without back reaction and showed that relativistic effects 
enhance mass accretion when back reaction is neglected. 

It is to be noted that one very important issue in under¬ 
standing the flow profile for accretion onto gravitating as¬ 
trophysical objects is the ‘transonicity’ of the flow. Let the 
instantaneous dynamical velocity and the local acoustic ve¬ 
locity of a compressible fluid moving along a space curve pa¬ 
rameterized by r be u{r) and o(r) respectively. Local Mach 
number M(r) of the fluid can then be defined as the ra¬ 
tio of the dynamical flow velocity to its sound speed, i.e.. 


M(r) = 


<^(r) 


a(r) 


. The flow will be locally subsonic or su¬ 


personic according to M(r) < 1 or > 1, i.e., according to 
u(r) < a(r) or u{r) > a(r). The flow is transonic if at 
any moment it crosses M = 1. This happens when sub¬ 
sonic to supersonic or supersonic to subsonic transition takes 
place either continuously or discontinuously. The points 
where such crossing continuously takes place are called sonic 
points and where such crossing takes place discontinuously 
are called shocks or discontinuities. One crucial difference 
between the flow characteristics around a black hole and 
around any other astrophysical object is, as the flow ap¬ 
proaches to any type of accretor other than a black hole, it 
can ‘physically’ hit the surface of the accretor directly and 
at the particular moment it collides with the surface, ac¬ 
cretion can be either subsonic or supersonic depending on 
the location of the sonic point (which itself is a function of 
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various accretion parameters) as well as on the location of 
the surface of the accretor whereas for accretion onto black 
holes, flow must fall onto the black hole only supersonically 
because even for the steepest possible equation of state, the 
maximum possible sound velocity attained will always be 
lower than the bulk velocity of the flow at the event hori¬ 
zon. On the otherhand, it is quite possible that the flow at a 
sufficiently large distance away from the accretor, would be 
subsonic in general (except for a special case when the su¬ 
personic wind from nearby object(s) falls onto the accretor). 
So black hole accretion is necessarily transonic to satisfy the 
inner boundary condition at the event horizon whereas ac¬ 
cretion onto other class of astrophysical objects may not 
always be transonic. 

One standard method to study the classical transonic Bondi 
flow is to formulate the basic conservation equations for the 
flow and then to simultaneously solve these conservation 
equations to get sonic quantities as function of various accre¬ 
tion parameters and also to calculate the values of various 
dynamical as well as thermodynamic quantities as functions 
of various accretion parameters or radial distance (measured 
from the central accretor in the unit of Schwarzschild ra¬ 
dius). However, the most effective approach to study the 
transonicity of spherically symmetric black hole accretion, 
as we believe, is to investigate the dependence of the loca¬ 
tion of the flow sonic point on various accretion parameters 
as well as to study the variation of Mach number of the flow 
with radial distance; which, perhaps, has not been explicitly 
discussed in existing literature for full general relativistic 
description of accretion onto black hole. In this paper, we 
would like to address this issue by formulating and solving 
the required equations of motion for a spherically symmetric 
transonic polytropic fluid accretion in Schwarzschild metric. 

2 GOVERNING EQUATIONS 

We take the Schwarzschild radius rg to be: 


2GMbh 



where Mbh is the mass of the black hole, G is universal grav¬ 
itational constant and c is velocity of light in vacuum. We 
assume that a Schwarzschild type black hole spherically ac¬ 
cretes fluid obeying polytropic equation of state. The density 
of the fluid is p(r), r being the radial distance measured in 
the unit of Schwarzschild radius rg. We also assume that the 
accretion rate is not a function of r and we ignore the self¬ 
gravity of the flow. For simplicity of calculation, we choose 
gravitational unit where unit of length is scaled in units of 
rg , unit of velocity is scaled in units of c and all other physi¬ 
cally relevant quantities can be normalized likewise. We also 
set G = c = 1 in system of units used here. 

For a Schwarzschild metric of the form 

ds^ = dt^ ~ ~ d.r^ “ D {d9^ -I- siv?6d(}P') 

the energy momentum tensor for a perfect fluid can be 
written as (Shapiro & Teukolsky 1983): 

rT-fOcQ a B , / a 3 oc0\ 

1 = eu u +p [u u — g ) 

where e and p are proper internal energy density and pres¬ 
sure of the fluid (evaluated in the local inertial rest frame of 


the fluid) respectively and m“ is the four velocity commonly 
known as 

^ ~ ds 

Equations of motion which are to be solved for our purpose 
are, 

1) Conservation of mass flux or baryon number conservation: 

(pMo). = 0 (la) 

and 

2) Conservation of momentum or energy flux (general rel¬ 
ativistic Euler equation obtained by taking the four diver¬ 
gence of T“^): 

{e + p)Ua- 0 U^ = -p,a - UaP,f 3 u'^ ( 1 &) 

where p is the proper matter density and the semicolons 
denote the covariant derivatives. 

Following Michel (1972), one can rewrite eq. 1(a) and eq. 
1 (b) for spherical accretion as 

Anpur^ = Min (2a) 

and 

as two fundamental conservation equations for time indepen¬ 
dent hydrodynamical flow of matter on to a Schwarzschild 
black hole without back-reaction of the flow on to the metric 
itself. Min being the mass accretion rate and C is some con¬ 
stant (related to the total enthalpy influx) to be evaluated 
for a specific equation of state. 

It is well known in general theory of relativity that a sta¬ 
tionary and axisymmetric space-time is endowed with one 
space-like [(^]) ] and one time-like [(^l) ] Killing field 
where (j> is the standard azimuthal co-ordinate. For a sim¬ 
pler space-time which is spherically symmetric, only the time 
like Killing filed is of particular interest and corresponding 
to this held, one can obtain the integral of motion along a 
streamline as: 

£ = heb^ ?-^eb (3) 

P 

where £, h (= and tb are the conserved specihc energy 
of the how including its rest mass, specihc enthalpy and 
specihc binding energy respectively. For polytropic equation 
of state p — Kp'^ (where 7 is the polytropic index and K 
is a constant which can be considered as the measure of 
the entropy of the how), the generalized expression for the 

sound velocity a = gives (Frank et. al. 1992, Wein¬ 

berg 1972): 



where T is the how temperature, p is the mean molecu¬ 
lar weight and niH^rnp is the mass of the Hydrogen atom. 
The subscript S indicates that differentiation is performed 
at constant specihc entropy. Also one can easily show that 
the specihc enthalpy h of the how is related to the sound 
speed through the following equation 
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Using the above expression for enthalpy, one can easily 
rewrite the conservation equation ( 2 b) as the specihc energy 
conservation as 


7 - 
-(l + a^) 



(6a) 


1 1 

Defining A4 = K to be another constant of 

motion for a shock free polytropic flow, one can rewrite eqn. 
(2a) as (see Chakrabarti (1996) and references therein): 


M 


dvr 


(7 - 1) 

7 — (a2 + 1) 




U€br 


( 66 ) 


The above equation may be considered as the outcome of 
the conservation of mass and entropy along the flow line. 
One can now easily derive the expression for velocity gradi¬ 
ent (by differentiating eq. 6 (a) and 6 (b)) as 

^ _ ujl- V?) [g^ (4r - 3) - l] 

dr 2r (r — 1 ) {v? — a^) 

Since the flow is assumed to be smooth everywhere, if the 
denominator of eq. 7(a) vanishes at any radial distance r, 
the numerator must also vanish there to maintain the conti¬ 
nuity of the flow. One therefore arrives at the so called ‘sonic 
point conditions’ by simultaneously making numerator and 
denominator of eq. 7(a) equal to zero and the sonic point 
conditions can be expressed as follows 


Uc — CLc — 




1 


4r,-3 


(76) 


here suffix c indicates that the values of the respective quan¬ 
tities are measured at the sonic point of the flow. For a 
specific value of £ and 7 , location of sonic point can be 
obtained by solving the following equation algebraically for 

Tc 


64r® - 1 ) - 16r^ “ 9) +4rc - 27) -7 27 = 0 

(7c) 

where ip = ( ) ■ It is important to note that though eqn. 

7(c) is a third order polynomial in r^, for all values of £ and 
7 , it gives only one real physical root for in general, we 
will return to this issue in next section. 

The spherical surface of radius r = can be defined as 
‘acoustic horizon’ because for r < r^, u > a\ and any 
acoustic disturbances created in this region are advected to¬ 
wards the black hole. Thus no acoustic disturbances created 
within this radius can cross the acoustic horizon and escape 
to the region r > r^- 

To determine the behaviour of the solution near the sonic 
point, one needs to evaluate the value of (^) at that point 
(we denote it by (^) ) by applying L ‘Hospitals’ rule on 
eq. 7(a). It is easy to show that (^) can be obtained by 
solving the following quadratic equation algebraically: 

/duO (7 - 1) (l 6 rc - 16rc - 87 rc-7 67 -7 3 ) rdu\ 
\dr), 3 r,( 4 r,- 3 )^ \dr). 
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Fig. 1 : Variation of the location of the accretion sonic point 
with conserved specific energy £ (which includes the rest mass 
energy of matter) and polytropic index 7 of the flow. While is 
plotted along z axis, £ and 7 are plotted along x and y axis respec¬ 
tively. The curves shown in the figure are plotted corresponding to 
the increments of 0.005 and 0.025 in £ and 7 respectively. Note 
that rc is produced even for flow with 7 < 4/3 and 7 > 5/3, 
though not all sonic points shown here allow real physical tran¬ 
sonic accretion. See text for detail. 

It is now quite straightforward to simultaneously solve eq. 
6 (a) and eq. 6 (b) to get the integral curves of the flow (curves 
showing the variation of Mach number with radial distance) 
for a fixed value of £ and 7 . Detail methodology for this 
purpose will be discussed in §3. 


3 SOLUTION PROCEDURE AND RESULTS 

One can obtain the location of the flow sonic point by solving 
eqn. 7(c) for a fixed value of the conserved specific energy 
£ and polytropic index 7 of the flow. As already mentioned, 
only one physical value of would be obtained which lies 
outside the event horizon. We solve eqn 7(c) for a range of 
values of £ and 7 for which a real physical solution is pos¬ 
sible. There are cases where only one real solution for is 
present but it lies inside the event horizon so we exclude 
those solutions. In Fig. 1., we show the variation of (plot¬ 
ted along z axis) with £ (plotted along x axis) and 7 (plotted 
along y axis). It is observed that the location of the sonic 
point non-linearly anti-correlates both £ and 7 which 
implies that ultra-relativistic flow u with low total specific 
energy will produce the sonic point located furthest distance 
away from the black hole. One thing is very interesting in 


(7 “ 1) (2)"c — 1) (24rc — 28rc — 8rc7 -7 4rc7 -7 87 -7 6 ) ^ 


2rc (4rc-3)5 (dr^ - 3) (r^ - 1) 


{Id) 


* Hereafter, we will describe the flow to be ultra-relativistic for 
7 = ^ and purely non-relativistic for 7 = ^ according to standard 
practice (Frank, et. al 1992). 
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Fig. 2: Integral curves of motion for a fixed value of ^ = 1.01 and 
7 = 4/3/. While Mach number of the flow M is plotted along y 
axis, distance from the central accretor (in units of Vg) is plotted 
along X axis in logarithmic scale. 

Fig. 1. We observe that sonic point is produced even for flows 
with 7 < 4/3 and 7 > 5/3 though not all such sonic points 
allow steady physical transonic flows passing through them, 
which would be clear from the following discussion and from 
Fig. 3. 

As we have mentioned above, Vc can be obtained from eqn. 
7(c). Now the value of (^) at Tc can be obtained by solv¬ 
ing eqn. 7(d) for a particular value of £ and 7 . Fourth or¬ 
der Runge-Kutta method is then employed to integrate from 
sonic point to get the dynamical flow velocity and the acous¬ 
tic velocity (so also the Mach number of the flow) at any 
point of the flow. It is well known that for spherically sym¬ 
metric accretion in Newtonian potential, two solutions are 
obtained while solving the conservation equations, one out of 
which corresponds to the accretion and the other one to the 
‘self-wind’. Same kind of situation is obtained here also for 
flows in general relativity. In Fig. 2 we plot variation of Mach 
number (plotted along y axis) with radial distance from the 
accretor (along x axis in unit of rg) in logarithmic scale for 
accretion (ABC) and wind (DBF) branches for a fixed value 
of £ (=1.01) and 7 (=4/3). The location of the sonic point 
B comes out to be 39.53 rg. It is important to note that 
though from eqn 7(c), sonic points are obtained for a wide 
range of values of £ and 7 , in reality, not all sonic points al¬ 
low a real physical flow to pass through them. We find that 
for accretion passing through some of the sonic points (ob¬ 
tained for a specific region of parameter space spanned by 
£ and 7 ), dynamical velocity or acoustic velocity, or both, 
becomes superluminal hence causality relation is violated 
which indicates that only a subset of the sonic points shown 
in Fig. 1 can be considered to study the transonicity of a 
real physical accretion flow. Here we would like to empha¬ 
size that still some sonic points, which allow real physical 
flow passing through them, are obtained for accretion with 
7 < I or 7 > |, see Fig. 3. 


Fig. 3: Behaviour of Mach number of the flow M (solid lines) 
and the flow temperature T (dotted lines) for a set of values of 
£ and 7 marked in the figure. M®’’’ and T®’' for the values of 
M and T respectively measured at a distance O.Olrg away from 

the event horizon. T®’’’ is scaled as T®’' -> 0.22 /T\\ (where 

Til = 10^^ °K) to fit in the same graph; see text for details. 

We would now like to investigate the behaviour of flow vari¬ 
ables close to the inner boundary of accretion. For this pur¬ 
pose we would like to study the the Mach number of the flow 
and the flow temperature very close to the event horizon of 
the black hole. As all equations diverge on the event horizon, 
we would like to tackle the problem in the following way: 
Let Th be the actual location of the event horizon (in units 
used here, rn = rg) and let Ar be equal to Srg where d is a 
small number less than one. We then define r®‘® = r;i-|-Ar to 
be the extreme point for our calculation of any flow variable 
Vfiow so that V/ioiu refers to the value of Vfiom measured at 
a radial distance (1 -|- <5) rg. Here we would like to calculate 
M®’” and T®’® for S = 0.01 for a fixed values of £ and 7 and 
to study the variation of M®”’ and T®’® with £ and 7 . Lower 
value of 5 can also be taken but it would only increase the 
computational cost, the general solution profile would re¬ 
main unaffected. 

In what follows, for a particular value of £ and 7 , we solve 
eqn 6 (a) and 6 (b) simultaneously upto r®“’ to get M®’® and 
in the same way we calculate M®‘® for various values of £ 
and 7 which allows real physical solutions. Similarly, as the 
sound speed a can be calculated easily while calculating M, 
one can obtain T®“’ (using eqn. (4)) for various values of £ 
and 7 . In Fig. 3 we show the variation of M®’” (solid line) 
and T®“’ (dotted line) with £ for a set of equi-spaced val¬ 
ues of 7 shown in the figure. T'”” has been normalized as 
T®^ —> 0.22 T®“’/Tii (where Tn = 10“ °K) to fit in the 
same graph. One can note that the value of minimum energy 
£m,in for which real physical accretion solution is available, 
non-linearly increases with increase of 7 . We see that M®’” 
non-linearly anti-correlates with £ and 7 . This is obvious 
because as we have seen from Fig. 1, flow with low energy 
and low 7 produces the sonic point far away from the black 
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hole so that the flow becomes more and more supersonic 
as it approaches the event horizon. However, follows 
exactly the opposite profile, it non-linearly correlates with 
both £ and 7 because for high energy accretion as well as 
for high enthalpy accretion (which is equivalent to high 7 
flow), flow temperature becomes higher. It is important to 
note that for a lower value of 5, both and keeps 
increasing keeping the general vs (£”, 7 ) profile 

remain unaltered. 


Weinberg, S. 1972, Gravitation and Cosmology : Principles and 
Applications of the General Theory of Relativity (John Wiley 
& Sons) 


4 CONCLUSION 

In this paper we have solved the equations governing general 
relativistic spherically symmetric hydrodynamic accretion of 
polytropic fluid on to a Schwarzschild black hole to investi¬ 
gate some of the transonic properties of the flow. We could 
study the variation of flow critical points with the conserved 
specific energy £ and the polytropic index 7 of the flow and 
showed that while high energy purely-non-relativistic accre¬ 
tion produces sonic points closer to the event horizon of 
the black hole, low energy ultra-relativistic flows pushes the 
sonic point far away from the hole. Also it is shown that not 
all sonic points allows a real physical flow to pass through 
them. Perhaps the most important finding in this paper re¬ 
garding the transonicity of the flow is that it is possible to 
obtain real physical sonic points (which allows a real phys¬ 
ical general relativistic spherically symmetric transonic ac¬ 
cretion through them) even for accretion with 7 < | or 

7 > |. Also we study the behaviour of the flow temper¬ 
ature and Mach number of the flow in close vicinity of the 
event horizon of the black hole as a function of £ and 7 and 
show that while low energy ultra-relativistic flow produces 
high Mach number at event horizon, high energy purely-non- 
relativistic flow produces high temperature at the horizon. 
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